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The synthesis of buckminsterfullerene [|[] has materiahzed the truncated icosahedral car- 
bon molecule, possessing the highest point group symmetry, and exhibiting pronounced 
Jahn- Teller (JT) effects 0] due to the electron-phonon coupling inherent in the high sym- 
metry. These effects determine many important physical and chemical properties of this 
third carbon isomer. Optical observations of phonon progressions and vibronic peak broad- 
ening in neutral Ceo in solution, matrix and solid films, as well as in its anion and cation 
counterparts clearly show the effect of JT distorted Ceo niolecule ^j^. On the other hand, 
the ESR lineshape also hints at the dynamic nature of these JT distortions 0,^. 

In order to study the JT distorted structures of Ceo ^ various theoretical ap- 

proaches have been used and quantum chemical CNDO/MINDO type calculations |§] 
have been performed. Based on group-theoretical and topological analysis, Ceulemans et al. 
have proposed that the symmetry of the stationary ground state of Cqq and CIq is de- 
termined by a weak 2nd-order coupling. Depending on the sign of the relevant parameter, 
the symmetry of this state may be Ds^ or Ds^. These two nearly degenerate configurations 
were also found in the quantum chemical calculations ||^. 

Meanwhile, within the Bogoliubov-de Gennes (BdeG) formalism, well-known in the liter- 
ature of conducting polymers [0, a simple SSH-like Hamiltonian is adopted to deal with the 
multi-mode JT problem of the molecule and the D^d symmetry structures of Ceo^ and Ceo" 
have been obtained |Tl|JT^ . For this static JT problem, the dynamic matrix of the different 
symmetry configurations has to be positive-definite to guarantee the dynamical stability. It 
has been found that the symmetry is not dynamically stable within the BdeG formalism 



T4 . The self-trapped excitons in neutral and charged Ceo niolecules have also been studied 
using this framework This self-consistent formalism is very efficient in describing the 
static JT distortions of this type of multi-electron, multi-mode systems. However, such a 
quasi-static approach is not sufficient. In fact, there are six equivalent D^d polarons obtained 
by breaking the Ih symmetry. These six polarons states are degenerate in energy and are 
connected by quantum tunneling, which may give rise to splitting of the energy state and 
restoring the full Ih symmetry. In fact, the binding energy of the D^d polaron is rather 



small (of the order of 40 meV as seen from Ref. [|T2| and the present calculation), which is 
comparable with the zero point energy for some of the JT active hg modes. Therefore, it 
is important to consider the dynamic JT effects due to quantum fluctuations. The genuine 
ground state should be an irreducible representation of the full Ih group and play a crucial 
role in related physical and chemical processes. 

In this paper, we will first briefly review the earlier results of the BdeG formalism, and 
describe the static solutions of the JT problem for Cgo" and Cqq anions. Then the first-order 
degenerate perturbation theory is adopted to determine the effect of the quantum fluctua- 
tions on the static configurations. The secular equations driving the tunneling splitting are 
solved analytically. We find that the quantum tunneling indeed leads to a further increase of 
the polaron binding energy by ~ 15 meV and the full Ih symmetry is recovered. Finally, the 
possible relations of our results with experiments and some neglected effects are discussed. 



In Ref. [|T2[, we simulated dynamic evolutions for a SSH-like system to find the stationary 
ground state configurations of Cgo and CIq molecules. We found that the singlet stationary 
configuration of Cqq also has D^d symmetry. To study the effect of quantum fluctuations 
on the static JT distorted D^^ polaronic configuration, we write down the full Hamiltonian 
in the second quantized form as 

H = He + -f^e-ph + -f^ph , (1) 
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where 



= E (f - C') • !7^(4.c,. + c]gc,g) . (2) 

<ij>,s 

Here a is the linear mode-independent vibronic coupling constant, is the amplitude of the 
normal mode in the homogeneous ground state as evaluated in Ref. |T^, is 



the unit vector along the unrenormahzed bond < i,j >. is the annihilation (creation) 

operator of the /xth normal phonon mode, defined as 6^|0 >= 0, in which |0 > is the un- 
renormahzed phonon vacuum. Each of the six D^a vacua breaks the Ih symmetry and can 
be defined by an origin-shifting operator C/^ as \o^ >— U^\o>, acting on the undistorted 
vacuum. This origin shifting unitary operator is constructed to translate the unrenormahzed 
phonon vacuum of lattice vibration to one of the six phonon vacua given by D^^ polaronic 
configurations {Q^}, P — 1, . . . , 6, 

U'-eM-li:]l^{K-W} ■ (3) 

Here Q^'s are c-numbers determined self-consistently through the gap equation — 
— < G/slF^lep >, where the multi-electron wavefunction je^j > is the solution of the 
corresponding BdeG equation in the same {Q^} configuration. 

The shifted phonon operator of /3th D^a polaron has the form of = and the 

corresponding phonon vacuum is defined as b^\oi3>— 0, where 

\o,>- exp(E ^KQ>> exp(- E ^iQif) (4) 

which satisfies the normalization condition, < o^|o;3 >= 1. b^^ describes the phonon exci- 
tation on the shifted vacuum. In the above description, we have already made an approx- 
imation, assuming that the D^^ polaron and the undimerized state have the same phonon 
spectrum, which will not significantly affect our result. 

The quantum tunneling between six minima is a tunneling process of the self-consistent 
lattice plus TT-clcctron states. Taking the states at the six equivalent D^^d polaronic minima of 
the adiabatic potential as the zero-order degenerate states, the state vector of the /9th polaron 
can be decomposed into two parts >— |e^ > |o/3 >, /3 = 1, 2, . . . , 6. The tunneling 
state should be a combination of the six non-orthogonal Dr^a polaron multi-electron states 
with different phonon vacua >— Z]^=iC^i|^/3 >, and can be solved by the degenerate 
perturbation theory. The secular equation is given by: 

||<*^/|i7|*^> <*^/|*^>|| = , (5) 
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where 



< V|if|*^> = < V|ife|^/3> + < Vl^e-ph|^/3> + < Vl^ph|^/3> , (6) 



<^^/|ife|^/3> = -to <ef}'\{cl,Cjs + c],Cis)\ef3><Of3>\0(3> , 

<i,j> s 

<Vl^e-ph|^/3> =a ^ <efs'\F^'\efs> hqf^ + <dp>\dp> , 

<i,j> s 



Here 



<0,,|5,>=exp{-X:^(Q?-O'} ' (7) 
and the matrix elements of < e/3/ 1 (cl^Cj s + c^j s'^is)\^p> and < C/?/ 1 e^j > are Slater determinants 



which can be calculated in accordance with the formulas listed in the Appendix of Ref. pi 
The overlap integral has the form of < \E'/3/|\E'/3 >=< e/3'|e/3 >< 0/3'|5/3 >. It is worth pointing 
out that <\E'^/|if|\E'/3>=<\E'/3|if|\l//3/>. 

To explicitly obtain the matrix elements, we number the axis pointing to the north pole as 
axis 1, while the remaining five axes spanning a pentagon-solid angle are numbered as 2 to 6, 
respectively. Due to symmetry consideration, we only need to calculate three kinds of matrix 
elements. When = /?, the diagonal matrix elements < \E'/3/|iJ|\E'/3' > (/3 = 1, 2, ■ ■ ■ , 6) are 
all equal, being the potential energy of D^^i polarons. For the Cgg" problem, the off-diagonal 
matrix elements only have two different values depending on the angle between the two 
axes. It is clear from the geometry that all the angles between two axes have a value ~ 63.4° 
except angles between axes 2 and 4, 2 and 5, 3 and 5, 3 and 6 and 4 and 6, which have a 
value ~ 180° — 63.4°. So it is only necessary to evaluate the elements < \Efi|if|\E'2 > and 
< \I'2|i/|\E'4 >. Moreover, since the HOMO state in D^^ polaron has an ungerade parity 
which is singly occupied, the Z^s^ polaron has an ungerade parity as well. Since the angles 
of these two pairs of axes are supplementary to each other, <\E'2|-fZ|^4>= — <^i|i/|^2>- 
The above analysis can also be applied to evaluate < \i//3/|\E'/3 >. The secular equation (^, 
now becomes 
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where a = Hn - E, t = U - ES, and Hu =< >, U =<^i|/J|\&2> and 
S =< \E'i|\[^2 >• Eq. (ID can be analytically reduced to 



{{Hn-Ef-5{U-ES)' 







(9) 



It is easy to obtain two three-fold degenerate eigenvalues, 



(10) 
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The six degenerate D^d polaron states span a six-dimensional space carrying a reducible 
representation of the Ih group, and the tunneling between them will remove this degen- 
eracy and reduce the six-dimensional space to two three-fold irreducible representations, 
T2u®Tiu, where the higher is identified to be levels by calculating characters of 
the representation matrices, while the lower E+ corresponds to Tiu levels. Moreover, the 
eigenvectors can also be determined analytically in terms of C/3i, which are listed in Table 
1, where the six tunneling split states |$j >, i = 1,2, ... ,6, are orthogonal to each other 
through the relations < $,1% >= 'E/s'isCpiC/sj < \E'/3'|\E'/3 >. From Eq. (0), we find the 
energy splitting, E^-E+ = (SHu - f/)2v^/(l - 55^). Analyzing this expression, we find 
that the factor SHu — U is actually equal to — < '^2\Hint\'^i >, where Hint = H — Hq 
and Hq is the zeroth-order Hamiltonian which has |\l'i> or |^2> as its lowest eigenstates. 
< \I'2|i^mi|^i > represents the hopping integral between the two states, which should be 
proportional to the energy splitting. Hence, we have 

V5 



E_ -E, 
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where —2 < \E'2|-ffmt|^i > is the energy sphtting for the two D^a polaron problem, and 
is an enhancement factor due to the presence of 5 neighbours. For the param- 
eter set of to = 1.91 eV,a = 5.0 eV/A, we got Hu = -206.66327 eV, U = -63.9725 eV, 
S = 0.30949, <^2|^mt|*i>= -0.01166 eV, = -206.579 eV, E+ = -206.679 eV, and 
an energy sphtting — Ej^ = 0.100 eV. 

We have also studied the Cqq polaron problem. The difference between a Cqq D^^ 
polaron state and a C\q D^d polaron state is that the former has a gerade parity, where the 
singlet HOMO electronic state is doubly occupied, whereas the latter has an ungerade parity. 
Therefore < \E'2|i/|\l'4 >=< ^i|if|^2 > for C^q polaron. The six-dimensional reducible 
representation can be reduced to Ag if^ as a result of tunneling splitting. The solutions 
of tunneling states and corresponding wavefunctions of C^q are given in Table 2. For the 
parameter set of to = 1-91 eV, a = 5.0 eV/A, we got Hu = -206.26237 eV, U = -9.0269 eV, 
and S = 0.043731, which resuhs in E„ = -206.255 eV, E+ = -206.290 eV, and an energy 
splitting, E^-E+ = 0.035 eV . 

According to the above results of the dynamic problem, the ground vibronic states of C^q 
and Cqq anions are ^Ti^ and ^Ag states, respectively, which are the same type of symmetries 
as the unrenormalized electronic states in the corresponding highest symmetric configura- 
tions of CIq and Cl^ . These results are consistent with the expected restoration of the 
initial 1^ symmetry through the quantum tunneling. In principle, the vibronic reduction ef- 
fect in these vibronic ground states is observable, especially in the ESR measurements. The 
estimate of Kato et al. ^ shows that the observed 0.003 reduction of g factor cannot be ex- 
plained without the inclusion of JT effect. However, the combined effects of JT vibronic and 
spin-orbit couplings need further theoretical and experimental investigations. Furthermore, 
the mixing of the vibronic ground states ^Ti„ and ^Ag with the low-lying excited '^T2u and 
^Hg states in the magnetic field should also be taken into account to determine the above 
observables. On the other hand, since the JT active modes hg are not dipole-active, the 
IR transition between the vibronic states ^Ti^ and ^T2„, ^Ag and ^Hg are parity forbidden. 
However, the tunneling splitting may be observable in the electronic resonant Raman scat- 
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tering measurements, if the degenerate low-lying excited ^T2„ and ^Hg states are taken as 
the intermediate states. 

Within our considerations, the pathways on the adiabatic potential surface from one D^a 
minimum to another have not been determined explicitly due to the complications result- 
ing from the multi-mode nature of the problem. Based on the work in the Ti„0(eg + t2g) 
problem |jl5|, Khlopin et al. ^ have shown that for the reduced linear one- mode Tiu<S>hg 



problem, the potential surface possesses a 2D equipotential continuum of minima in the 
strong coupling limit, and indicated that the dynamic properties cannot be solved by ne- 
glecting the multimode nature via a simple scale transformation because of the dynamic 
interaction between the scaled coordinates. On the other hand, the locations and heights of 
the possible D^^, and C2/1 saddle points on the potential surface, corresponding to the 
lower ranking epikernels of Ih-, are the key factors to set the possible pathways. However, 
unfortunately, the knowledge of the height of these saddle points is not sufficient to fully 
determine the pathways, because the quadratic vibronic couphng is absent in our present 
approach. 

Another neglected factor in the above discussions is the effect of electron-electron cor- 
relations, which could affect both static and dynamic aspects of the problem through the 
competition with the vibronic coupling fl^, especially, in the case of multiply-charged an- 
ions. However, the evidence of the important role of vibronic coupling in the the dynamic 
process of the optical excitation experiments is clear enough to show that the lattice of Cgo 
is dynamically distorted and the symmetry is well defined (not as -Dsd) The dynamical 

JT effects of charged fuUerenes have recently been discussed within a very different approach 
(Berry phase) [1^. The large tunneling splitting obtained in this calculation supports the 
assumption of spherical symmetry, adopted there. 

In summary, the static linear Ti^ Shg JT effect is analyzed for the C^q and Cqq an- 
ions based on the numerical results within the BdeG formalism. The dynamic problem has 
been solved analytically within the standard perturbation approach. The quantum tunnel- 
ing among the six stationary equivalent D^d polaronic minima results in the ^Ti„ and ^Ag 
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vibronic ground states, accompanied by the low- lying ^72u and ^Hg vibronic states for C^q 
and Cqq anions, respectively. Therefore, in experiments one should look for signatures of 
these states rather than the Dq^ states obtained by neglecting the dynamical JT effects. 
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Table 1, Energies and Wave Function Coefficients Cpi of Tunneling states in C, 



Wave function Coefficients 
Sym. E N{Cii, ^2%, C^i, C^, C^i, Cg^) 



(10(1 - y55))-i/2(y5, -1, -1, -1, -1, -1) 
< (40(1 - V55))-V2(o, 4, -1 - v^, -1 + v^, -1 + v^, -1 - V^) 
(4(5 - v^)(l - v^>5))-V2(o, 0, -1 + v^, -2, 2, 1 - V5) 



lu 



1+V55 



(10(1 + a/5>S))-V2(V5, 1, 1, 1, 1, 1) 

(40(1 + ^5))-i/2(0, 4,-1 + ^5,-1-75,-1-75,-1 + ^5) 
(4(5 + 75)(1 + 75^))-V2(o, 0, 1 + 75, 2, -2, -1 - 75) 
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Table 2, Energies and Wave Function Coefficients Cp of Tunneling states in C|( 



Wave function Coefficients 
Sym. E Ar(c,,C2,C3,C4,C5,C6) 



(30(l-5))-V2(5, -1, -1, -1, -1, -1) 
(20(l-5))-V2(o, 4, -1, -1, -1, -1) 
< (12(l-5))-V2(0, 0, 3, -1, -1, -1) 
(6(l-^))-V2(o, 0, 0, 2, -1, -1) 
(2(l-^))-V2(o, 0, 0, 0, 1, -1) 

v 



A, ^ (6(l + 5^))-V2(i, 1, 1, 1, 1, 1) 
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